Muon-spin-rotation study of magnetism in Na x Co02 single 

crystals with 0.78 < x < 0.97 

C. Bernhard 1,2 , Ch. Niedermayer 3 , A. Drew 1 , G. Khaliullin 2 , S. Bayrakci 2 , 
J. Strempfer 2 , R.K. Kremer 2 , D.P. Chen 2 , CT. Lin 2 , and B. Keimer 2 

1 ) University of Fribourg, Department of Physics 
and Fribourg Center for Nanomaterials (FriMat), 

Chemin du Musee 3, CH-1700 Fribourg, Switzerland. 

2) Max- Planck- Institute for Solid State Research, 
Heisenbergstrasse l,D-70569 Stuttgart, Germany, and 

3) Paul-Scherrer-Institut, CH-5232 Villigen PSI, Switzerland. 
(Dated: February 6, 2008) 

Abstract 

By muon spin rotation we investigated the magnetic properties of a series of highly Na doped 
Na-cCoC^ single crystals with 0.78(1)< x <0.97(1). Our data provide evidence for an intrinsically 
inhomogeneous magnetic state which can be described in terns of hole doping (Na vacancy) induced 
magnetic clusters that percolate at 1-x > 0.04 until they yield a bulk magnetic state near x =0.78. 
Evidence for a strong (likely geometrical) frustration of the magnetic order is obtained from the 
anomalous doping dependence of the spin fluctuation rate (above the ordering temperature) which 
is strongly enhanced at x =0.78 as compared to x =0.97. 
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The discovery of super conductivity with T c <5 K in the hydrated layered cobaltate 
Nao.3sCo02 * 1.3 H 2 |l| has renewed the interest in the physical properties of the par- 
ent compound Na x Co02- Its rich phase diagram includes several phases with extraordinary 
electromagnetic and thermoelectric properties . One attractive feature is the triangu- 

lar coordination of Co which may favor geometrical frustration and unconventional ground 
states js, 0|. Another curious aspect concerns the finding that the anomalous electromag- 
netic properties are not confined to the Mott-insulator regime at low x, but also occur in 
the vicinity of the anticipated band-insulator at x=l. Examples are the occurrence of mag- 

tinnnn □ 

netic order at x > 0.75 [3, |8L |9j, [10|, 111], signatures of strong spin-charge coupling [I2j and 



giant thermoelectric effects [4j . The nature of the underlying magnetic state is still debated. 
While initial data were interpreted in terms of an incommensurate spin- density- wave (SDW 
recent neutron- and muon-spin-rotation (/xSR) studies on single crystals 0, 9, 10, 11] 
yield commensurate A-type antiferromagnetic (AF) order and a weakly anisotropic magnetic 
exchange coupling constant that contrasts with the large Fermi-surface anisotropy. Besides, 
an interpretation in terms of localized Co moments in the low-spin-state is also questionable 
since the long range magnetic order would occur in a strongly diluted state with only 1-x 
magnetic Co 4+ ions ( 5 t 29 , S=l/2) and x non magnetic Co 3+ ions ( 6 t 2g , S=0). 

Here we present muon spin rotation (/xSR) and dc magnetization measurements which 
provide new insight into the evolution of the magnetic properties of highly Na doped 
Na x Co02 single crystals with 0.78(1)< x <0.97(1). Our data provide evidence for hole 
doping (or Na vacancy) induced magnetic clusters which remain isolated at 1-x <0.04 but 
percolate at l-x>0.04 until a homogeneous magnetic state is achieved near 1-x ~0.25. We 
also find that this long-range ordered magnetic state exhibits signatures of strong (likely 
geometrical) frustration. 

The Na x Co02 crystals have been cleaved from ingots that were grown with an optical 
floating-zone furnace [jjj]. The average Na:Co ratio of our crystals has been carefully de- 
termined by inductively coupled plasma spectroscopy (ICPS) and by neutron activation 
analysis with an accuracy of about 1%. The crystals were also characterized by electrical 
transport, dc magnetometry, and single crystal x-ray diffraction. The latter suggest that the 
x=0.78 crystals are single phase (likely at), while the ones at x=0.87, 0.92 and 0.97 contain 
mixed a- and a- phases, similar as was previously reported 0, 4]. For the latter, the inten- 
sity of the at- phase peaks was found to exhibit a strong variation, even for crystals with a 
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similar bulk Na content (according ICPS and neutron activation analysis). Accordingly, we 
suspect as outlined below that the x-ray data my be strongly affected by Na deficient sur- 
face layers. Finally, the /iSR measurements were performed with the GPS setup at the 7rm3 
beamline of the Paul-Scherrer-Institut (PSI) in Villigen, Switzerland, which provides 100% 
spin-polarized muons. Details concerning the technique can be found in [15j. Here we only 
note that zero-field /iSR technique provides information about the local scale distribution of 
the internal magnetic fields of magnetic materials. In particular, it yields a reliable estimate 
of the bulk volume fractions of the magnetic phases, even for the case of weak and strongly 
disordered magnets. 

First we discuss the //SR data of the Na .97CoO 2 crystal with l-x=0.03(l) holes per Co 
(Figure 1). The zero- field (ZF) spectra between 5 and 300 K are shown in Fig. 1(a). The 
low-T spectra evidently consist of two distinct parts. Nearly 40% of the signal depolarize 
rapidly while the remaining 60% exhibit a much slower relaxation. This suggests that 
the sample is in a spatially inhomogeneous magnetic state with about 40% of the volume 
containing sizeable magnetic moments and around 60% remaining non magnetic. It is also 
evident that the correlations in the magnetic regions persist to temperatures in excess of 
200 K. The absence of oscillations in the fast relaxing signal implies a broad distribution 
of the magnetic fields that originates either from static but spatially disordered moments 
or else from dynamical fluctuations. It is well known that a large longitudinal field (LF) 
enables one to distinguish between these two cases 15]. It gives rise to a strong reduction 
of the amplitude of the magnetic signal (so-called "decoupling effect") in the static case but 
not in the dynamic one. Accordingly, the data at LF=1 kOe in Figs, lb and lc establish 
that a transition from a static to a dynamic state takes place between 15 and 30 K. The 
signatures of the corresponding spin freezing transition are also apparent in the transverse 
field (TF) relaxation rate, A TF , in Fig. Id. The slow increase of X TF between 300 and 
50 K is consistent with a gradual decrease of the spin fluctuation rate, while the following 
rise and the saturation below 20 K are characteristic of a freezing transition at Tf ^20 K. 
Furthermore, the transition width and the large value of Xq F are characteristic for a glassy 



state 16] . This freezing transition is also evident in the dc magnetic susceptibility in Fig. 2 
due to a deviation from the Curie law below 50 K and a cusp around Ty=20 K, as marked 
by the arrow. 

The inhomogeneous magnetic state for samples with x ~ 1 was previously noticed but 
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it was discussed in terms of chemical inhomogeniety due to inclusions of a Na deficient and 
thus magnetic phase [itJ or even of impurity phases [3]. The former interpretation was 
motivated by x-ray diffraction data which signal the presence of two Na x Co02 phases with 



different Na content 



0, fl 



181 ] . However, in the following we argue that such a purely 



chemical scenario is not sufficient to explain our data. Instead, we present evidence that the 
formation of intrinsic, doping induced magnetic nanoclusters seems to play an important 
role. For example, we investigated crystals from two growth batches with x =0.97(1) out of 
which one also exhibited the x-ray signatures of a secondary Na deficient a-phase, while the 
other according to x-rays was an almost pure a-phase (see inset of Fig. Id). Despite of this 
substantial difference, our /iSR data yield virtually identical bulk magnetic properties for 
both crystals (Fig 1 shows our /iSR data of the pure a-phase crystals). The only distinction 
concerns a minor difference in the magnetic volume fraction that is well consistent with the 
1% uncertainty in the Na content. These observations suggest that one needs to distinguish 
between a nanoscale phase separation as probed by /zSR and a macroscopic one as evidenced 
for example with x-rays. While a strong correlation between the Na vacancies and the doped 
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es which determine the electronic and magnetic properties of the C0O2 layers is evident 
191, it remains unknown how the nanoscale magnetic clusters are related to the macroscopic 
Na deficient regions that a r e fre^ent.v observed in s-ray flQ and „ent r o„ diffraction 
We noticed however, that the x-ray signature of the a-phase in highly Na doped crytstals 
appears to be time dependent, even on the scale of weeks. To the contrary, we confirmed that 
the /xSR data do not exhibit any noticeable changes, not even after one year. Accordingly, 
we suspect that the apparent instability of the crystal surface at ambient conditions towards 
the formation of Na 2 C03 and a subsequent development of a Na deficient surface layer may 
play an important role here. We note that this is especially relevant for polycrystalline 
samples on which the majority of experiments has been performed llll, llSj. In this 

context, we emphasize that our crystals are mm sized and that /xSR is a truly bulk sensitive 
technique since the muon implantation depth is of the order of 100 /im. 

Our combined /iSR and dc magnetisation measurements are indeed consistent with the 
point of view that the relevant length scale of the magnetic phase separation at x=0.97 is on 
the order of nanometers. In the first place they are suggestive of an antiferromagnetic spin 
coupling within the magnetic clusters. Specifically, the dc magnetisation data in Fig. 2 yield 
a fairly small value of the total Curie-moment is otherwise difficult to reconcile with the large 
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magnetic volume fraction and the sizeable magnetic Co moment that emerges from the /iSR 
data. The solid lines in Fig. 2 show fits to the high temperature magnetisation data for T > 
50 K with the function, Xmoi = jSq+X°- The Curie- Weiss term with C = (l—x)NAi~tljf/3k B , 
Avogadro number, , effective moment, // e // , and Boltzman constant, kg, accounts for 
the paramagnetic component. The T independent term, x°, accounts, besides Van-Vleck 
paramagnetism and core diamagnetism, for the strong correlation effects within the clusters 



as detailed in 



201 ] . The parameters obtained at x=0.97 for the field directions parallel (||) and 



perpendicular (_L) to the c-axis are 0|| i j_=-61 and -45 K, xjj ± = 1-10 - 4 and 0.78-10" 4 cm 3 /mol 
and C|| i j_=8.4-10 -3 and 5.5-10 -3 cm 3 /mol, respectively. For 3% doping this translates into a 
moment per doped hole of // e // = 1.5 and 1.2 in units of Bohr magneton, fiB, respectively. 
With a spin only Lande factor of g = 2 we thus derive a total spin of Sim = 0.4 and 0.28 
per cluster. Such a small magnetic moment per cluster needs to be reconciled with the 
sizeable moment per magnetic Co ion of about 0.1-0.3 \i& that has been inferred from /xSR 
and neutron experiments Bay , While this value was obtained for AF samples with 



x ~ 0.82, our /iSR data indicate that the Co moments at 2=0.97 are of similar magnitude. 
The local field at the muon site, B M , as deduced from the TF-/iSR depolarization rate of 
\q F ~ 23 at 2=0.97 (assuming randomly oriented moments) of <B fM >=^— « 270 G 
(7 M = 851.4 MHz/T is the muon gyromagnetic ratio) compares indeed well to the one 
obtained from the highest precession frequency of the antiferromagnetic samples of =3.3 
MHz (see Fig. 3) with Bn—^^ ~ 240 G. It is now interesting to note that the non 
magnetic regions experience noticeable magnetic stray fields which apparently originate from 
the magnetic clusters since they also undergo a freezing transition around 20 K. This is 
evident from the ZF- and LF-/zSR data in Figs lb and lc which show that the slowly 
relaxing component exhibits a decoupling effect for LF=50 Oe at 15 K that is absent at 30 
K. The large magnitude of the decoupling effect at 15 K suggests that stray fields essentially 
persist throughout the entire non-magnetic volume of the sample. This result needs to be 
compared with the estimates based on dipolar field calculations which indicate that the 
stray fields around antiferromagnetically ordered clusters are extremely short ranged, i.e. 
they are noticeable only over a few nanometers. Accordingly, our combined /zSR and dc 
magnetisation data are indeed suggestive of a nanoscopic coexistence of the magnetic and 
the non magnetic regions. 

This conclusion is supported by our finding that with increasing hole doping, 1-2, the 
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magnetic clusters percolate and form extended magnetic patches until a homogeneous mag- 
netic state is achieved near x=0.78. This is shown in Figure 3 which displays the ZF-/iSR 
time spectra for x=0.92, 0.87 and 0.78. The insets show the distribution in frequency space 
as obtained with a maximum entropy analysis. It is evident that all the 5 K spectra contain 
an oscillatory component due to magnetic order at least in parts of the volume. We used 
the function P(t) = P(0) Yl ^-i cos(27r^i) exp(— X l t) to obtain the parameters for the 

Fig. 4 together with previous data for £=0.82(1) [8j. Starting from a threshold of x tr ~0.04, 
the amplitude of the oscillatory signal increases rapidly with hole doping, and approaches 
100 % around £=0.78. The full oscillatory amplitude of the 5K spectrum at x=0.78 also 
confirms that the magnetic field at the muon site, B ^ is almost perpendicular to the muon 
spin direction (which is tilted by 60 degree against the c-axis). Accordingly, the /iSR am- 
plitudes should indeed be proportional to the corresponding magnetic volume fractions |{|. 
The c-axis orientation of B M furthermore suggests that the muons reside at interstitial sites 
near the oxygen ions rather than on vacant Na sites where B M would have a larger com- 
ponent parallel to the C0O2 layers and be much smaller in magnitude as was previously 
noted js, 11]. While we do not attempt here a detailed assignment of the muon sites, we 
notice that the signal at z/ M =1.2 MHz likely arises from the muons that are stopped near the 
boundaries of the magnetic patches, since it has a large relaxation rate of about 3 /xs _1 , it 
dominates at x=0.92 and 0.87 and is absent at x=0.78. Correspondingly, we suggest that 
the signals near 2.5 and 3.0 MHz, which become very narrow and further split at x=0.78, 
arise from muons that reside within the order patches. Whether the distinction between 
these muon sites arises from different structural (for example due to the Na neighbors) or 
magnetic environments requires further investigation. 

Next, we discuss the remarkable trend that the /xSR spin fluctuation rate, r M , at T>Tjv 
is strongly enhanced for x=0.78 as compared to x=0.97. At x=0.97 slow spin fluctuations 
with 10~ 6 s > T/j, >10 -10 s persist to at least 200 K while at x=0.78 the spin fluctuations 
become too fast to be detected within the /xSR time window of r /t <10 -10 s immediately 
above T N —22 K. This is evident from Fig. 3c where already the 25 K spectrum exhibits 
a Gaussian shape and a small relaxation rate that is entirely due to the nuclear moments. 
This trend is supported by the magnetization data for £=0.78 (Fig. 2) which exhibit a 
broad maximum around 60 K that is characteristic for low-dimensional and/or geometrically 
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frustrated systems with strong (quantum)fluctuations. We emphasize that this behavior once 
more is not compatible with a purely chemical phase separation where the spin fluctuations 
should be enhanced due to finite size effects for the isolated nanoscopic magnetic clusters 
at a;=0.97. Our data therefore highlight that the ordered magnetic state at x=0.78 is not 
a conventional kind of A-type antiferromagnet that should not be subject to such strong 
frustration effect (not even for a triangular lattice). 

In summary, by /xSR and dc magnetisation measurements we have explored the magnetic 
properties of a series of highly Na doped Na x Co0 2 single crystals with 0.78(1)< 2 <0.97(1). 
Our data provide evidence for an intrinsically inhomogeneous magnetic state with hole 
doping (Na vacancy) induced magnetic clusters that percolate at 1-2 > 0.04 until they yield 
a bulk magnetic state near x =0.78. We observed a distince magnetic provide evidence that 
the 

The SST may be induced by the doped holes which lower the crystal field symmetry of 
the neighboring Co 3+ ions and thus give rise to a splitting of the e g levels that reduces the 
energy difference to the highest occupied t 2s level. Accordingly, at low doping each hole 
induces a magnetic cluster with a central low spin (LS, S=l/2) Co 4+ ion that is surrounded 
by six intermediate spin (IS, S=l) Co 3+ ions 12|, 113J. The SST model predicts indeed that 
a homogeneous long-range ordered magnetic state occurs at 1-2=0.25 Q, Q|. Explicit 
calculations have shown that while the dominant interactions within the clusters are AF, 
a weak residual FM in-plane interaction between the clusters gives rise to an A-type AF 



state 
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20j . Notably the SST model predicts that the induced IS Co 3+ ions reside on 

nn 

a 2-d Kagome lattice [121 , 1131 ] which is geometrically frustrated and thus subject to strong 
quantum fluctuations. Furthermore, the SST model accounts for the rapid disappearance 
of magnetism at l-2>0.25 where the symmetry reduction associated with the doped holes 
becomes insufficient to induce a SST. It also accounts for the evolution of the magnetic 
volume fraction at low doping, although we note that this can be equally well explained 
by chemical phase separation into non magnetic NaiCo02 and magnetic Nao.85Co02 
At 1-2=0.03 the SST model predicts that each hole gives rise to a magnetic cluster which 



contains 24 oxygen ions and subsequent muon sites [8|, [11] that are adjacent to at least one 
magnetic Co ion. Accordingly, due to the 0:Co ratio of 2:1, it yields a fraction of 36% 
magnetic muon sites which agrees well with our experimental value of about 40%. The 
SST model also allows one to understand the glassy behavior and the strongly reduced spin 
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fluctuation rate of the isolated clusters in terms of a delicate energy balance between single-, 
double- or even triple hole clusters that likely depends on the spatial arrangement of the Na 
vacancies. The related disorder and the distribution in the shape anisotropy of the clusters 
thus may well account for the glassy magnetic properties at l-a;=0.03. 

We conclude by pointing out that a corresponding doping induced SST is indeed well 
established for the simple perovskite Lai-^Sr^CoOg whose Co ions are in a correspond- 



ing valence state and crystallographic coordination [2l| . Nevertheless, the possible SST in 
Na x Co2 has some unique aspects like the small total moment of the magnetic clusters of 
only fi e ff ~ 1.5/ie and the circumstance that the IS Co moments reside on a Kagome lat- 
tice and thus are geometrically frustrated. Both of them favor strong quantum fluctuations 
which may well be relevant for the unconventional magnetic, electronic and thermoelectric 
properties of these fascinating materials. 
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FIG. 1: The /uSR data for Nao.97Co02- (a) Displays the zero-field spectra at representative tem- 
peratures. The solid lines show fits as obtained with the function P(t)/P(0) = 0.4 • exp (—Ait/ + 
0.6 • KG ■ exp (— A2i) where KG represents the so-called Kubo-Toyabe function which accounts for 
the nuclear magnetic moments, (b) and (c) show longitudinal-field spectra at 30 K and 15 K, 
and (d) the transverse field (TF) relaxation rate of the fast relaxing component. Inset: X-ray 
diffraction data (Cux« with wavelength, A =1.54 A) for crystals from two different growth batches 
with x=0.97(l) indicating the presence of a pure a-phase or mixed a- and a' phases, respectively. 
Arrows mark the position of the a'-phase (002) peak and the a-phase (003) peaks. 
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FIG. 2: T-dependent magnetic susceptibility of Na^Co-O^ with x=0.97 and 0.78. Solid lines show 
Curie- Weiss fits for x=0.97 as specified in the text. The arrow marks the freezing transition at 
T f =20.7 K. 



FIG. 3: Zero-field juSR spectra for Na I Co02 crystals with x=0.92, 0.87 and 0.78 which exhibit a 
rapidly increasing fraction of the oscillatory signal below T/v- 
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FIG. 4: Doping dependence of the fit parameters of the magnetic component in the ZF-/iSR spectra 
at 5 K. (a) Amplitude of the precession frequencies. The half down (up) filled symbols show the 
amplitudes of the splitted frequencies at x=0.78 with 2.2 (2.6) MHz and 3.1 (3.3) MHz, respectively. 
Solid (open) circles show the sum of the oscillatory (total magnetic) amplitudes. The dotted line 
is a guide to the eye (b) Doping dependence of the corresponding relaxation rates as shown by the 
same symbols. 
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